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ABSTRACT 

This paper is the second in a series in which we present intermediate-resolution, wide- 
wavelength coverage spectra for a complete sample of 14 compact radio sources, taken 
with the aim of investigating the impact of the nuclear activity on the cirumnuclear 
interst ellar medium (IS M) in the early stages of radio source evolution. In the first 
paper (jHolt et al.l l2008f) we presented the kinematic results from the nuclear emis- 
sion line modelling and reported fast outflows in the circumnuclear gas. In this paper, 
we use the line fluxes to investigate the physical conditions and dominant ionisation 
mechanisms of the emission line gas. We find evidence for large electron densities and 
high reddening in the nuclear regions, particularly in the broader, blueshifted com- 
ponents. These results are consistent with the idea that the young, recently triggered 
radio sources still reside in dense and dusty cocoons deposited by the recent activity 
triggering event (merger/interaction). In addition, we find that the quiescent nuclear 
and extended narrow components are consistent with AGN photoionisation, split be- 
tween simple-slab AGN photoionisation and mixed-medium photoionisation models. 
For the nuclear broader and shifted components the results are less clear. Whilst there 
are suggestions that the broader components may be closer to shock plus precursor 
models on the diagnostic diagrams, and that the electron temperatures and densi- 
ties are high, we are unable to unambiguously distinguish the dominant ionisation 
mechanism using the optical emission line ratios. This is surprising given the strong 
evidence for jet-cloud interactions (broad emission lines, large outflow velocities and 
strong radio-optical alignments), which favours the idea that the warm gas has been 
accelerated in shocks driven by the radio lobes expanding through a dense cocoon of 
gas deposited during the triggering event. 

Key words: galaxies: ISM - Galaxies, galaxies: active - Galaxies, ISM: jets and 
outflows - Interstellar Medium (ISM), Nebulae, galaxies: kinematics and dynamics - 
Galaxies 



1 INTRODUCTION 

It is now understood that AGN-feedback, in the form of out- 
flows in the circumnuclear ISM, is likely to be a crucial part 
of galaxy evolution and is responsible for the close relation- 
ship betwe en the black hole mass and the galaxy bulge prop- 
erties (e.g. iFerrarese fe Merrittll2000l: iGebhardt et al.l l 2000l: 
iTremaine et al.ll2002l : iMarconi fe HuntfeOOot ). However, due 
to a lack of observational evidence, feedback is often in- 
voked as a 'black box' and often relies on powerful AGN 
winds to drive the outflows and shed the natal cocoon de- 
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posited during the triggering event (e.g. m erger / interaction; 
|Pi Matteo et alj2005UHopkins et al]|2005l ). Whilst this may 
be true for radio-quiet AGN, it is likely that in radio- loud 
AGN, the expanding radio jets also cont ribute to the feed- 
back through jet-induced outflows (e.g. IBicknell fe Dopital 
119971 ). particularly during the early stages of the radio source 
evolution when the radio source and ISM are on similar 
scales. 

Recent work (by e.g. iBest et alJkood : iBower et aill2006l 
fe lCroton et al . 2006) has shown that extended radio sources 
have a significant impact on the hot gas on cluster scales 
through feedback effects (e.g. the quenching of cooling flows 
in massive haloes). Hence, a significant contribution to 
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the overall AGN feedback from the young, expanding ra- 
dio sources on smaller scales should be expected. Indeed, 
over the last few years, fast nuclear outflows (up to ~ 
2000 km s" 1 ) have been detected in a number of compact 
(young) radio sources, both in the optical emission l ines (e.g. 
ITadhunter et al.ll200ll ; fHolt et al.|[2003l . 120061 . [20081 : the lat- 
ter 3 references: H03, H06 & H08 hereafter) and in neutral 
hydrogen as deep absorption features (e.g. iMorganti et al.l 
120031 , 12004 120051 ). 

In H08, we reported results for the emission line kine- 
matics in a complete sample of 14 compact radio sources, in- 
cluding both Compact Steep Spectrum radio sources (CSS: 
size < 15kpc) and Gigahertz-Peaked Spectrum radio sources 
(GPS: size < 1 kpc). The observed line widths and blueshifts 
(up to ~ 2000 km s _1 ) in the nuclear emission line gas 
are too large to be explained by gravitational motions. In- 
stead, the extreme kinematics were interpreted as signatures 
of strong interactions between the young, small-scale radio 
jets and the emissio n line gas, following the model proposed 
for PKS 1549-79 by iTadhunter et alj pOOlh . This interpre- 
taion is supported by evidence for strong radio-optical align- 
ments in all CSS sources with de ep HST imaging, inc l uding 
some sources in this sample (e . g. de Vries et al.| [l997. 1999; 
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Although the expanding radio jets provide a convenient 
outflow driving mechanism, consistent with the results on 
the emission line kinematics, it is important to test this sce- 
nario further. In this paper we use optical emission line ratios 
(from the line modelling in H08) to investigate the phys- 
ical conditions and the dominant ionisation mechanism(s) 
of the narrow line gas. In particular to determine whether 
the emission line regions are shocked, as expected in the 
case that the o utflows are being accelerated by jet-cloud in- 
teractions (e.g. IVillar-Martfn et al.lll99Sl : I Clark et all 1 19981 : 
ISolorzano-Inarrea et al.ll200ll ). In Section 2 we give a brief 
overview of the sample, observations, data reduction and 
analysis techniques and refer the reader to H08 for a detailed 
discussion. In Section 3 we discuss the line fluxes and, in par- 
ticular, the reddening, density and temperature of the emis- 
sion line gas, giving brief details of the line modelling. We 
discuss the line ratios in detail in Section 4, along with diag- 
nostic diagrams and radio-optical correlation plots. All line 
fluxes, aperture information and extracted nuclear spectra 
are presented in the Appendices. Information on the emis- 
sion line kinematics for all sources is presented in H08. 

Throughout this paper we assume the following cosmol- 
ogy: H = 71 km s~\ Q = 0.27 and Q A = 0.73. 



2 SAMPLE SELECTION, OBSERVATIONS, 
DATA REDUCTION AND ANALYSIS 
TECHNIQUES 

This paper is the second in a series of two reporting deep 
optical spectroscopic observations of a complete sample of 
compact radio sources. Below we summarise the important 
points and refer readers to H08 for a detailed description of 
the sample, observations, data reduction and analysis tech- 
niques. 

Our complete sample comprises 14 compact radio 
sources derived from the northern 3C and 4C samples and 
southern 2Jy sample and includes 8 CSS, 3 GPS, 2 compact 



core sources and 1 compact flat-spectrum source. The sam- 
ple has intermediate redshifts (z < 0.7) and a radio power 
range of 26 < log P 5G h z < 28 (W Hz" 1 ). 

Long-slit optical spectroscopic observations of the full 
sample were obtained during three observing runs with 
ISIS/WHT, EMMI/NTT and FORS2/VLT. In order to in- 
clude the outflowing regions in the slit, and to ensure the 
spectra were of sufficient resolution, the observations were 
made using a 1.0-1.5 arcsec slit, while to reduce the effects 
of differential atmospheric refraction, all objects were ob- 
served either at low airmass (sec z < 1.1) and/or with the 
slit aligned along the parallactic angle. Due to various ob- 
servational constraints, we have only aligned the slit along 
the radio axis for approximately half of the sources. For full 
details of the observational setup, we refer readers to Table 
2 of H08. 

The data were reduced in the usual way (bias subtrac- 
tion, flat fielding, cosmic ray removal, wavelength calibra- 
tion, flux calibration) using the standard packages in iraf. 
To ensure good flux calibration (±5 per cent), several spec- 
trophotometric standard stars were observed during each 
run with a wide (5 arcsec) slit. This accuracy was confirmed 
by good matching in the flux between the red and blue spec- 
tra. □ Additional standard stars were observed with a narrow 
slit to correct the spectra for atmospheric absorption fea- 
tures (e.g. A and B bands). The spectra of all sources were 
corre cted for Galac ti c ext inction usin g the E(B-V ) values 
from ISchlegel et ail (|l99Sf ) and the ISeatonl (|l979l ) extinc- 
tion law. The wavelength calibration accuracies were 0.06- 
0.53A, 0.06-0.24A and 0.20-0.24A for the WHT, NTT and 
VLT data respectively, dependent on wavelength range. Sim- 
ilarly, the spectral resolutions were typically 3.3-4.8±0.2A, 
4.3-6.7±0.lA and 6.5-7.4A for the WHT, NTT and VLT 
data, again dependent on wavelength range. The spatial res- 
olutions were 0.36 arcsec/pixel, 0.33 arcsec/pixel and 0.25 
arcsec/pixel for the WHT, NTT and VLT data respectively. 
The seeing range was 0.8-2.5 arcsec. Full details are given in 
Table 2 of H08. 

As discussed in H08, in order to detect and accurately 
model the broader emission line components, it is necessary 
to model and subtract the underlying continuum emission. 
For the broader, blueshifted emission line components, vari- 
ations in the continuum can make it difficult to establish the 
centroid, FWHM and line flux of the Gaussian. If the un- 
derlying continuum also contains strong stellar absorption 
lines (e.g. H/3 etc), a good model of the continuum is re- 
quired to ensur e that the emission line fluxes are accurately 
modelled. iHoltl (120051 ) estimated that subtracting different 
continuum models which all provided a reasonable fit to the 
SED can cause a difference of up to 10% in the flux in the 
broader components, although typically the narrower com- 
ponents are less affected. 

As a first step in modelling the continuum, we gen- 
erated and s u btrac ted the nebular continuum following 
iDickson et alj (|l995l ). taking full account of reddening fol- 
lowing the techniques outlined in H03. The remaining 
continuum was then modelled using a customised idl 
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data used is for lines observed with the same arm of the ISIS 
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minimum \ 2 fitting programme ( see e .g. iRobinsonl l200ll; 
iTadhunter et all 12003 ; iHolt et all 120071 for details) allow- 
ing up to three separate continuum components: 12.5 Gyr 
unreddened Old Stellar Population (OSP), a Young Stel- 
lar Population (YSP) with reddenings < E(B-V) < 1.6 
and age 0.01 < £ysp < 5.0 Gyr, and an AGN/power law 
component. The best fitting model, defined as that with the 
least number of components required to adequately model 
both the overall SED and discrete stellar absorption features 
(e.g. Ca H+K, Balmer series), was then subtracted. The 
continua of three sources were not modelled in this way: 3C 
277.1 (quasar), PKS 1814-63 (foreground star) and 3C 303.1 
(mis-matching of the blue and red spectra). Full details of 
the subtracted models can be found in H08. 

The spectra were extracted and analysed using the 
STARLINK packages FIGARO and DIPSO. Note, some compo- 
nents for a few key diagnostic lines were not detected in 
some apertures. For these, we have estimated upper limits 
for the flux by adding Gaussians to the data until the line 
component could just be detected on visual inspection of the 
spectra. 

Throughout this paper we use the kinematic component 
definition defined in H08: 

• narrow: FWHM < 600 km s -1 ; 

• intermediate: 600 < FWHM < 1400 km s _1 ; 

• broad: 1400 < FWHM < 2000 km s" 1 ; 

• very broad: FWHM > 2000 km s _1 . 



3 RESULTS 

As discussed above, the nuclear emission lines in compact 
radio sources are often broad, with asymmetric profiles re- 
quiring multiple Gaussian components to model them. In 
H08, we modelled the nuclear emission lines for all sources 
in our sample, identifying between 2 and 4 Gaussian com- 
ponents with velocity widths and blueshifts of up to 2000 
km s _1 . Here, we use the emission line models from H08 
and focus our analysis on the emission line fluxes and ra- 
tios. We use emission line ratios and diagnostic diagrams to 
investigate whether the compact radio sources in our sample 
still retain their dense and dusty natal cocoons (reddening 
and emission line density), and to investiage the dominant 
ionisation mechanism(s) and the implications for the likely 
outflow driving mechanism. For brevity, the line fluxes of all 
lines in all apertures extracted are presented in Appendix 
A. 

3.1 Aperture extraction and emission line 
modelling 

For a detailed discussion of the emission line modelling, we 
refer readers to H08 and we only summarise the key points 
here. 

3.1.1 Aperture selection 

In H08, with the exception of modelling the most extended 
emission lines in the spatial direction to accurately deter- 
mine the systemic velocities, the main results on the emis- 
sion line kinematics were derived from wide nuclear aper- 



tures centred on the nuclear continuum emission. These nu- 
clear apertures, which also form the basis of the analysis 
in this paper, are shown in Appendix [B] along with spatial 
profiles of the bright emission lines and the continuum emis- 
sion, and the full extracted nuclear spectra. In addition, 5 
sources (3C 268.3, 3C 277.1, PKS 1345+12, PKS 1814-63 & 
3C 459) show evidence for extended emission in several lines 
and we have extracted further apertures. The positions of 
these, along with the extracted spectra, are also shown in 
Appendix B. 

3.1.2 Emission line modelling 

In general, we have modelled all emission lines in a partic- 
ular aperture using a kinematic model derived from the [O 
III]AA4959,5007 doublet, a standard technique for modelling 
the emission lines in radio galaxies. In short, the velocity 
widths and shifts of the various components were fixed to be 
the same as for the [O III] doublet, with further constraints 
arising from atomic physics. As all lines can be reproduced 
by the same model, we can assume with confidence that the 
emission originates from similar regions of the ISM. Hence, 
it is meaningful to use emission line ratios to investigate the 
physical and ionisation properties of the emission line gas. 

Notable exceptions are the nuclear apertures of the GPS 
sources PKS 1345+12 and PKS 1934-63. For these aper- 
tures, one model did not reproduce all lines. In these latter 
cases, only components consistent in all lines (e.g. the nar- 
row components in PKS 1345+12) are plotted in the diag- 
nostic diagrams. Finally, for 3C 303.1, we only consider line 
ratios from lines in the same arm of the instrument as the 
two arms do not match in flux. 



3.2 Emission line densities 

During the early stages of radio source evolution, the nuclear 
regions of radio sources are likely to harbour a rich ISM, 
deposited during the activity triggering event. Here, we in- 
vestigate the density of the circumnuclear gas using the [S 
II] 6716/6731 density diagnostic ratio. The measured ratios 
were converted to densities using the iraf program temden, 
part of the Stsdas package, which is based on the five -level 
atom calculator developed bv lde Robertis et al.l (1987). For 
all calculations, we have assumed an electron temperature of 
T c = 10,000K, although the [S II] ratio varies little with tem- 
perature. As discussed above, when modelling the emission 
lines in the nuclear aperture of each source, the profiles of 
all lines were well fitted by the corresponding [O III] model. 
The only exceptions were PKS 1345+12 (discussed in detail 
in H03) and PKS 1934-63. 

The derived electron densities are summarised in Table 
[1] In the nuclear apertures, we have been able to accurately 
determine the electron density for 5 components only: 3C 
213.1 (narrow), 3C 303.1 (intermediate), PKS 0023-26 (nar- 
row), PKS 1934-63 (narrow) and 3C 459 (narrow). For the 
remaining components, we have only been able to deter- 
mine whether the components are consistent with low or 
high density. The reason for this is that the blend is highly 
complex, requiring 4 to 8 Gaussian components (2-4 per line) 
for an adequate fit, and the velocity widths and shifts of the 
broader components are often comparable to the separation 
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Figure 1. The [S II] AA6716,6731 doublet in the nuclear aper- 
ture of PKS 1345+12 (taken from H03). The faint solid line (his- 
togram) represents the data, the bold line is the overall best fitting 
model to the doublet, the dashed line is the [O III] model to the 
[S II] doublet and the dot-dashed, solid and dotted pairs of lines 
are the narrow, intermediate and broad components. This plot 
highlights the difficulties in obtaining a unique model for the [S 
II] doublet. Although it is relatively straightforward to fit the cor- 
rect number of components, and we can confidently determine the 
total flux in, for example, the intermediate or broad components, 
the widths and shifts of the broader components are comparable 
to the separation of the lines in the doublet. This makes it diffi- 
cult to determine the ratio between the two intermediate or two 
broad components, which is critical for determining the density. 

of the lines in the doublet (see Figure [1] and also H03) ; it 
is usually possible to confidently determine the total flux 
in e.g. the broad components together, but to determine the 
density, the ratio is required. Hence, whilst the [O III] model 
generally provides a good fit to the [S II] doublet, in many 
cases it is not possible to obtain a physically viable good fit 
without forcing the 6716/6731 ratio to the high (0.44) or 
low (1.42) density limit of the diagnostic. 

Although we have not been able to accurately measure 
the electron densities in the majority of apertures, from the 
results presented in Table [1] we can draw the following gen- 
eral conclusions: 

(i) The nuclear narrow components tend to be consistent 
with low electron densities (few hundred cm ), as do the 
extended narrow components. 

(ii) The nuclear broader components (intermedi- 
ate/broad/very broad) appear to be consistent with higher 
densities (n e > 10 3 cm -3 ). 



3.3 Reddening 

In Section f3. 21 and H03, we showed that there is a sugges- 
tion that the nuclear regions of our sample of compact radio 
sources may harbour a dense ISM (n e > few 1000 cm ). 
Observationally, in addition to high densities, large amounts 
of gas and dust can also be detected as reddening of the op- 
tical spectrum. 

We have investigated reddening in all apertures in 
our sample using the Balmer line ratios (e.g. Ha/H/3, 
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Figure 2. The Ha,[N II]AA6548,6584 blend in the nuclear aper- 
ture of PKS 1549-79. The faint solid line (histogram) is the data, 
the bold line is the overall best fitting model to the doublet and 
the dotted and dashed lines are the narrow and intermediate com- 
ponents (3 of each: 1 from Ho and 2 for [N II]). The solid Gaus- 
sian is the broad component detected only in Ha in the optical 
(see H06 for details). This plot highlights the complexity of the 
Ha/[N II] blend and the difficulties in measuring the reddening 
using Hcr/H/3 ratio. 

H7/H/3 and H5/H/3), assum ing a simple foreground screen 
model for interstellar dust (|Seatonlll979l ). The correspond- 
ing E(B-V) values for each object /aperture are summarised 
in Table □ 

Note, as discussed in Section 2 and H08, many galaxies 
in the sample exhibit stellar absorption features. Whilst the 
continuum modelling and subtraction attempts to remove 
these features, if a weak Balmer emission line is coincident 
with a strong absorption feature (e.g. H5 and, to a lesser 
extent, H7), this may result in an inc orrect (large) estimate 
of the reddening l|Fosburv et alJlioO ll ) . When modelling and 
subtracting the continuum (see H08), the effect of the un- 
derlying absorption lines on the emission lines was investi- 
gated, focussing in particular on the different components 
of H/3. The well defined narrow components of the emission 
lines are often relatively unaffected by the subtraction of 
different models. However, the flux in the broader compo- 
nents could be affected by as much as 10% or more. Hence, 
wherever possible, estimates of reddening from the stronger 
Balmer lines are used. It should also be noted that in the 
more extreme sources, with several components to the emis- 
sion lines, there is also difficulty in accurately modelling Ha, 
due to the complexity of the blend with the [N II] doublet. 
Figure [2] demonstrates this for the nuclear aperture of PKS 
1549-79. 

The degree of estimated reddening varies over a wide 
range, from negligible through to highly extinguished. There 
does not appear to be any particular trend with radio source 
size (i.e. GPS or CSS) although some of the highest E(B- 
V) values are measured in the GPS sources (e.g. E(B-V) = 
1.0 in 4C 32.44 and E(B-V) = 1.4 in PKS 1345+12). The 
main trend appears to be the increase in reddening with the 
FWHM of the emission line component within a particular 
object. In 7/14 sources (denoted by the stars in Table[T]), the 
degree of reddening increases significantly from the narrow 



Table 1: Summary of the physical conditions (density, reddening & temperature) for the sample. Where there are two components with similar 
FWHM (i.e. two narrow components), the component at the systemic velocity is presented first (for PKS 1306-09 the blue narrow component is 
presented first), t was measured and found to be consistent with either the low density limit ([S IIJA6716/A6731 = 1.42) or the high density limit 
([S IIJA6716/A6731 = 0.44), within the uncertainties. $ forced the the low (1.42) or high (0.44) density limit. See text for details. * sources in which 
there is a clear trend of increasing density (column (g)) or reddening (column (I)) with FWHM. Blank positions in the table indicate kinematic 
components not observed in a source/aperture. For components which arc not detected in a relevant diagnostic line, the entry is -. 
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(d) 


(e) (f) 
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(h) 


(i) 


(i) 


(k) 


(I) 


3C 213.1 


nuc 










zero 


>0.2 






* 


3C 268.3 


nuc 
SE 
NW 


1.42$/1.42$ 
1.42f 

_ 


1.42$ 






zero / zero 

zero 

0.618J 


0.7+81 

— 0.1 






* 


3C 277.1 


nuc 

SE blob 
SE HII 


<600 

-/- 


0.44$ 






zero 

zero / zero 
zero 


0-718:2 




BLR:zero 


* 


4C 32.44 


nuc 


-1- 




- - 




zero/zero 




0.9181 




* 


PKS 1345+12 


nuc 


1.42$ (<150) 


> 5300 


> 4200 


* 


zero 


o.4i8;i 


I 4+0-5 




* 




NW,PA160 


500128° 








zero 








SE,PA160 


1.42t/1.42$ 








o.5i8;i/i.2i8;^ 












NE,PA230 


1.42$ 








o.7l8:l 












SW,PA230 


-/1.42$ 








-/o.7i8:i 










3C 303.1 


nuc 


1.42$ 


iin+40 
1±u -40 




*? 


zero 


5+ ' 2 
U.O_ 2 








PKS 0023-26 


nuc 


4ooi?™; 


1.42$ 






zero 


zero 








PKS 0252-71 


nuc 








zero / zero 


zero 








PKS 1306-09 


nuc 










zero / zero 










PKS 1549-79 


nuc 


1.42$ 


0.44$ 






zero 


zero 




BLR:2.8l8;l 






PA75 


1.42$ 


0.44$ 






0.418:1 


4+ 01 








PA-5 


1.42$ 


















PKS 1814-63 


nuc 
ext 


1.42$/700l 1 1 ^ 
240+ 100 








-/- 










PKS 1934-63 


nuc 


600+ 300 




0.44$ 




zero 










PKS 2135-20 


nuc 




"/- 








zero/1.318'2 








3C 459 


nuc 
ext 


300+ 100 
1.42$ 




0.44$ 
1.42$ 




zero 
zero 


i - 3 -0.1 
zero 




* 
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Table 1. continued. 



Object Aperture Electron temperature, 10' 5 T e (K) 

n i b vb 

(a) (6) (I) (m) (n) (o) 



3C 213.1 


nuc 


<19 


<30 


3C 268.3 


nuc 
SE 
NW 


16i 8 3 /15^ 
18±\ 


18±i 8 


3C 277.1 


nuc 

SE blob 
SE HII 


25±\ 

<35/<55 

<95 




4C 32.44 


nuc 


11^/<18 




PKS 1345+12 


nuc 

NW.PA160 
SE,PA160 
NE,PA230 
SW,PA230 


u±l 




3C 303.1 


nuc 


<16 




PKS 0023-26 


nuc 


22+ 3 


<35 


PKS 0252-71 


nuc 


32t?°/30i?° 


110 +ioo 


PKS 1306-09 


nuc 


«W±S8/40±*g 




PKS 1549-79 


nuc 

PA75 

PA-5 


14+ 4 




PKS 1814-63 


nuc 
ext 






PKS 1934-63 


nuc 






PKS 2135-20 


nuc 




<18 


3C 459 


nuc 
ext 


<46 





component, often suffering negligible reddening (although 
some objects exhibit narrow components with moderate red- 
dening) , through to the intermediate and broad components 
where reddening is almost always significant. This effect was 
first observed in the GPS source PKS 1345+f2 and is dis- 
cussed by H03. Of the remaining 7 objects which do not 
show this effect, it is not possible to do this test in 2 sources 
(i.e. both components observed have similar FWHM as in 
PKS 1306-09) and in 1 further object, PKS 1814-63, the de- 
gree of reddening is not estimated. Taking this into account, 
7/11 sources (or 67 %) show the trend of increasing redden- 
ing with increasing FWHM. Hence, all reddening estimates 
are consistent with dense and dusty cocoons surrounding the 
young radio sources, and in more than half of the sources ob- 
served there is strong evidence for stratification of the ISM 
as in PKS 1345+12. Such strong reddening in the nuclear 
regions of radio sources is not uncommon and is observed 
in a number of extended radio sources s uch as Cygnus A 
l|Tadhunter et al.lll994l ; lTavlor et al.ll2003h . 

Whilst we can see the general trends in the reddening 
and draw conclusions about the existence of a dense and 
dusty circumnuclear cocoon, given the various uncertainties 
discussed above, we have not corrected the line fluxes for red- 
dening. For some line ratios (e.g. [S II] for the density), the 
lines are close together and reddening effects will be negli- 
gible. For the diagnostic diagrams presented later, for ratios 
in which the lines have a large separation in wavelength, we 
plot the uncorrected values. 



3.4 Electron temperature 

The electron temperature can be an important diagnostic for 
distinguishing between the different ionisation mechanisms. 
For example, one would expect to observe higher tempera- 
tures if the gas is shock ionised rather than photoionised by 
the AGN. 

The electron temperature, T e , can be derived from the 
temperature diagnostic [O III] (4959 + 5007) /4363. Where 
the weaker [O III]A4363 line is detected, we have calcu- 
lated the electron temperature, assuming an electron density 
of 10 3 cm -2 . The derived electron temperatures are sum- 
marised in Table [1] 

No clear trends are observed across the entire sample 
but all measured temperatures are relatively high (in the 
main > 14,000 K), and in some objects (e.g. 4C 32.44) a gra- 
dient with FWHM is observed, although the uncertainties 
are large. Such temperatures are too high to be accounted 
for by simple single slab photoionisation models which gen- 
erally predict temperatures of < 11,000 K (jTadhunter et all 
1989) and may be a signature of shock ionisation, although 
mixed-medium ioni sation models could also account for such 
high temperatures (|Binette et al.|[l996h . Density effects may 
also be responsible for the unrealistically high temperatures 
measured in some sources/apertures. 
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4 DISCUSSION 

4.1 Do compact radio sources harbour a dense 
and dusty natal coccoon? 

iTadhunter et all l|200ll 'l concluded that the compact, flat 
spectrum radio source PKS 1549-79 is a young, recently 
triggered radio-loud AGN which is still surrounded by a 
dense and dusty natal coccoon, obscuring the active nu- 
cleus from view. This was confirmed by the detection of 
the proto-quasar, in itially at near infra-red wavelengths 
( Bell amy et al.l[2003t ) and, later, at the reddest optical wave- 
lengths in higher quality optical spectra (H06) Deep optical 
images also reveal that PKS 1549-79 has recently been in- 
volved in a major merger event. Similarly, H03 reported the 
discovery of significant reddening and high electron densities 
in the nuclear aperture of the GPS source PKS 1345+12. 
Again these results were interpreted as the signature of a 
dense and dusty natal cocoon, particularly as PKS 1345+12 
also shows evidence of a recent merger (double nucleus in a 
common envelope e.g. lAxon et al.|[2000h . 

In Sections 13.21 and 13.31 we investigated the density 
and reddening in the nuclear apertures of 14 compact radio 
sources, including PKS 1345+12 and PKS 1549-79. Whilst 
it has been difficult to obtain unique results, particularly 
for the electron densities, we can draw some general con- 
clusions. The majority of the sources in the sample do tend 
to show evidence for high densities and large reddening in 
the nuclear regions, particularly in the broader shifted com- 
ponents. This trend is particularly apparent in the redden- 
ing - 7/14 sources show convincing evidence for increasing 
reddening with FWHM (2 GPS, 3 CSS & 2 compact core 
sources) . 

Hence, it is clear that the nuclear regions of compact 
radio sources contain large amounts of gas and dust. This is 
consistent with the larg e reddening observed in the optical 
spectra of CSS quasars (|Baker fc Hunst ead 1995), although 
at X-ray wavelengths, the evidence for differences in col- 
umn densities between y oung AGN and more ex tended radio 
sources is contradictory. iGuainazzi et al.1 (120061 ) report high 
column densities in GPS sources, for all sources in their small 
sample with good quality hard x-ray data. The measured 
column densities in GPS sources are higher compared to FR 
I and BLR FR II sources, but comparable to high-ionization 
FR II sources, and so they suggest the obsc uring mate- 
rial is located in an obscuring torus. However, IVink et al.l 
(2006) find a variety of column densities in GPS and CSS 
sources, with no evidence for higher column densities com- 
pared to more extended radio sources. This latter observa- 
tion is consistent with the fact that strong optical reddening 
is not uncom mon in the nuclear regions of extended radio 
sources (e.g. IRobinson et al.1 Il987l; Tadhunter et al.l Il994l ; 



iRobinson et al.ll2000l ; lRobinsonll200ll ; iTavlor et al.ll2003f )" 

Whilst there is evidence for a dense and dusty circum- 
nuclear cocoon, there appears to be insufficient material to 
constrain and fru strate the radio sources (e.g. PKS 1345+12: 
iHolt et alJliool ). The only way frustration would be possi- 
ble is if the matter were distributed i n giant clouds, whic h 
may be able to disrupt a radio jet (|Bicknell et al.l [2003) . 
Hence, the observations are consistent with compact radio 
sources being young, recently triggered radio sources resid- 
ing in dense and dusty natal cocoons, rather than old, con- 
fined and frustrated radio sources. 



4.2 What is the dominant ionisation mechanism 
for the gas? 

ITadhunter et al.1 (2001) suggested that the outflow in 
the compact radio source PKS 1549-79 was driven by the ex- 
panding radio jets. Certainly, jet-driven outflows are a con- 
venient way of explaining the extreme kinematics observed 
in the optical emission line gas. However, it is also possi- 
ble that the outflows have been accelerated by some other 
mechanism, for example, quasar-induced winds or starburst 
driven superwinds. 

A key way to try to distinguish between these scenarios 
is the use of emission line ratios and diagnostic diagrams, a 
technique successfully employed for identifying jet-cloud in- 
teractions in the ext ended emission line r egion s of extended 
radio sources (e.g. IVillar-Martm et all 1 19981 ; IClark et al.l 
Il998l : Is"olorzano-Inarreall200ll '). We have therefore plotted a 
variety of diagnostic diagrams for the different emission line 
components in Figure0 (nuclear narrow & extended compo- 
nents) and Figure [4] (nuclear broader/shifted components). 
In addition to the data for this sample, we have also plotted 
a comparison sample of nearby extended radio galaxies with 
detailed line ratio information (see the caption of Figure 3 
for details) and model results for simple slab AGN p hotoioni- 
sation (from MAPPINGS Ic: seelFerruit et al.|ll997h . mixed- 
medium photoionisation (|Binette et al.lll996T ) and the lat- 
est sh ock and shock plus precursor models from lAllen et al.l 
(2008). Full details of the models are given in the caption 
of Figure [3] Note that it can be difficult to distinguish be- 
tween the AGN and shock plus precursor models for the high 
shock velocity case, because the hard continuum generated 
in the hot post-shock gas (that photoionizes the precursor) 
can mimic that of an AGN nucl e us. M oreover, the mixed 
medium models of iBinette et al.l (|l996h , while successfully 
explaining the relatively strong high ionization lines and 
high electron temperatures measured in some AGN, cover 
only a limited range of parameter space for the matter- and 
radiation-bounded components included in the models; bet- 
ter agreement with the observed line ratios in a given object 
might be obtained by tuning the input parameters of such 
models. 

Initial inspection of Figures [3] and [4] reveals that the 
positions of the broader and narrower emission line compo- 
nents of the compact radio sources are not significantly dif- 
ferent from each other, and are remarkably similar to those 
of the comparison sample of extended radio sources. This 
latter result is surprising given the striking differences in the 
nuclear emission line kinematics and the strong radio-optical 
alignments reported for CSS sources (see Section 4.3). 

In terms of the results for the compact radio sources, 
whilst the diagrams are complex, it is still possible to iden- 
tify some tendencies in the data when interpreting all dia- 
grams together: 

• Narrow components. In the majority of the diagrams 
(e.g. Figures 3a, 3b, 3e), the nuclear narrow components lie 
close to tracks for photoionisation (both simple AGN pho- 
toionsition and the A M /j sequence) whilst in Figures 3c & 
3d, the nuclear narrow components lie away from the plot- 
ted photoionisation tracks, and appear more consistent with 
shock models. However, in the latter cases, the data always 
lie closer to the higher velocity end of the tracks (v Bhock > 
500km s _1 ). Given that the nuclear narrow components are 
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(a) 

Figure 3. Diagnostic diagrams for the entire sample of nuclear narrow components (large red open triangles) and extended components 
(large green open triangles). Arrows represent upper limits. The ionisation models plotted are as follows: 

— crosses joined by solid black lines: optically thick power-law (AGN) photoionisation sequence in the ionisation parameter U (2.5 X 
10~ 3 < U < 10 _1 ) from mappings for a=-1.0, -1.5 & -2.0; 

— bold black solid line: mixed-medium models (including both matter-bounded and ionisation-bounded clouds; 10~ 2 < 5: 10) from 
iBinette et all [I1996T) ; 

— dotted lines: new shock models from [Allen et all (120081 ). Sequences are for magnetic paramter, B/y'n = 0.01, 1, 10, 100, 1000 fiG 
cm 3 / 2 , with shock velocities 10 2 — 10 4 km s _1 and pre-shocked gas density n e = 100 cm -3 . The stars identify the highest shock velocity 
(10 4 km s — 1 ) on each track with the large open star identifying the tra ck with the highest magnetic paramter; 

— dashed lines: new shock plus precursor models from I Allen et all ||2008|) . Sequences are for magnetic paramter, B/y^n = 0.01, 1, 10, 100, 
1000 /iG cm 3 / 2 with shock velocities 10 2 -10 4 km s — 1 and pre-shocked gas density n e = 100 cm -3 . The stars identify the highest shock 

velocity (10 4 km s -1 ) on each track with the large open star identifying the track with the highest magnetic paramter. 

For comparison, various data from the literature ar e overplotted (small b l ack points): nuclear regions (filled triangles: IWills et al] 
20021; IStorchi-Bergmann et all Tl99rj; iTadhunterl Il987l; IRobinson et all Il987l, iGrandi fc Osterbrockl |l97St , ICostero fc Osterbrockl Il97jl; 
Cohen fc Osterbrockl Il98ll : iDicksonl Il997l . ISolorzano-Inarreal l200ll: IVillar-Marti'n et al.lll998l: IClaxk et*aDll998l: iRobinson et al.l|200(fl : 
extended emission li ne regions (EELRs, filled squares: iTadhunter et all Il994l ; iMorgan ti et al l 1 199 it IStor chi-Bergmann et al. 1996: 
Robinson_et^ al. 2000); EELRs with evidence for jet-cloud interactions (open squares: Villar-Martm et all dl998h : IClark et all dl998f) : 
S^lOTzano-Karreal ll200lh ) and all regions in Cygnus A (larger black open symbols: Taylor et al. 2003). 



relatively quiescent, with small velocity widths (FWHM < 
400km s _1 ), and no blueshifts (H08 demonstrated that they 
are at the systemic velocity), such an interpretation would 
be unlikely. It should also be noted that, in both Figures 3c 
& 3d, the locus of the data for the narrow components is 
similar to that for the mixed-medium photoionisation mod- 



els - varying the parameters of the mixed-medium mod- 
els would result in tracks which describe the data better 
than the other models. We therefore interpret Figure 3 as 
the nuclear narrow components showing a tendency to be 
consistent with AGN photoionisation (both simple slab and 
mixed-medium models). The extended narrow components 
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are also plotted on Figure 3, and we see no significant differ- 
ences between the nuclear and extended narrow components 
suggesting they have similar ionisation mechanisms. 

• Broad, shifted components. As for the nuclear nar- 
row components, careful consideration of all information to- 
gether, along with details of the models, is required in order 
to disentangle the overlapping models on the diagnostic di- 
agrams in Figure [4] From all the diagrams in Figure 4, we 
can rule out pure shock ionisation. In all 5 panels of Figure 
4, the data occupy the same region as the shock plus precur- 
sor models, typically in the region of high velocity (v s hock > 
500 km s _1 ) shocks (this is most clearly seen in e.g. Figures 
4a, 4b & 4c). In some diagrams, the location can also be 
interpreted as consistency with AGN photoionisation (e.g. 
Figures 4a & 4b), but in Figures 4c & 4d, the data is nei- 
ther in the same region as, nor has the same locus as, any 
of the photoionisation models. Again, by combining all dia- 
grams, there appears to be a mild tendency towards shock 
plus precursor models, with high shock velocities (v B hock ^ 
500 km s _1 ) in the nuclear shifted components, although 
this evidence is not conclusive. 

• Object classification. The [O III] A5007/H/3 vs. [N 
II]A6584/Hq diagram (Figures 3e & 4e for the narrow and 
broad components respectively) is a well known diagram to 
distinguish between the different classes of obj ect and there- 
fore the dominant ionisation mechanism (e.g. iKewlev et al.l 
l200ll ; Kauffma nn et all 12003). Almost all components are 



consistent with AGNs/Seyferts, with a few in the region of 
the diagram occupied by LINERS. There are a couple of 
notable exceptions - the extended region to the SE of the 
nucleus in 3C 277.1, and the extended regions to the SE 
and SW in PKS 1345+12 are all consistent with stellar pho- 
toionisation. The observation of ongoing star formation in 
radio galaxies is rare, with only a few other cases known: 
e.g. Com a A (ISolorzano-Inarrea fc Tadhunterl [20031 ) ; PKS 
1932-46 (|Villar-Martfn et al.ll2005h . The Super Star Clus- 
ters (SSCs) in the halo of PK S 1345+12 are discussed in 
iRodriguez Zaurin et "all (|2007t ) . 



4.3 What drives the outflows? 

In H08 we reported evidence for extreme emission line out- 
flows in the nuclear regions of 11/14 compact radio sources. 
From the kinematics alone, the large velocity shifts and 
widths (FWHM) observed (up to 2000 km s _1 ) are en- 
tirely consistent with the idea that the young, expanding 
radio source is strongly interacting with the natal cocoon 
and driving outflows in the ISM, as suggeste d by the re- 
sults for PKS 1549-79 (|Tadhunter et al-lfeoOll ). This argu- 
ment is further strengthened by the growing evidence in the 
literature that the line emission in CSS sources is both on 
similar scales to, and often closely aligned with, the radio 
source (|de Vries et al.lll997l . [l99st I Axon et~al]|2000l ). This 
effect is observed for all CSS sources with HST imaging. Re- 
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Figure 5. Radio-optical correlation plots. Plotted are the sources from this sample: CSS (red triangles), GPS (blue squares), compact 
flat spectrum (magenta diamonds) and compact core (cyan diamond) . The comp arison sample is a redshift selected subsample (0.05 < 
z < 0.7) of the the 2 Jy sample (e.g. iTadhunter et al.lll993l ; iMorgant i et al,| [l993) and is plotted as small filled triangles (FR II sources), 
small filled squares (FR I sources) and crosses (core-jet sources). The plots are as follows: (a) [O III]5007 luminosity versus 5GHz radio 
power; (b) [O IIJAA3727 luminosity versus 5GHz radio power; (c) H/3 luminosity versus 5GHz radio power; (d) [O II]/[0 III] ratio versus 
5GHz radio power and (e) [O III]/H/3 ratio versus 5GHz radio power. 



cent HST imaging of two of the sources in our sample ( PKS 
1345+12 and PKS 1549-79) bv lBatcheldor et al l (|2007r ) has 
also shown that in the smaller GPS sources, the optical line 
emission is also on a similar scale to the radio source. How- 
ever, although galaxy wide outflows (i.e. starburst-driven 
winds) can be confidently ruled out for GPS sources, on the 
basis of emission line region morphology it is not possible to 
distinguish between AGN- winds and jet-cloud interactions, 



which would both be expected on the nuclear scalefl Whilst 
the kinematics and radio-optical alignments are clearly sug- 
gestive of jet-cloud interactions and radio jet-driven out- 
flows, the above evidence alone is, so far, circumstantial. 

In this paper, we have attempted to address this prob- 
lem using optical emission line ratios and diagnostic dia- 

2 Due to the small spatial scales, it is currently not possible to 
resolve the outflowing regions sufficiently. 
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grams. In Sections 3.2 and 3.4, we measured electron densi- 
ties and temperatures. Although we measure large densities 
and high temperatures (n e up to ~ few 1000 cm -3 '; T e > 
14,000K) in the broader, blueshifted components, both of 
which are expected for shocked gas, due to the large mea- 
surement uncertainties, we cannot say with certainty that 
the densities and temperatures in the broader components 
are significantly different from those in the narrower com- 
ponents. 

In Section 4.2 and Figures [3] and |4j we plotted a variety 
of diagnostic diagr ams, including the latest shock ionisa- 
tion models from lAllen et all (|2008T ) along with some of the 
older AGN photoionisation models (fro m mappings and th e 
mixed-medium (A M /j) sequence from iBinette et ail [l996) . 
However, as discussed above, the evidence from the diagnos- 
tic diagrams is far from clear, with only a weak trend for the 
broader components to be more consistent with fast (v s h oc fc 
> 500km s _1 ) shock plus precursor models. 

Our results are therefore rather surprising. The evidence 
in the literature for jet-cloud interactions in the nuclear re- 
gions of compact radio sources, based on the emission line 
alignments and extreme emission line kinematics, is strong. 
However, at best, we have found only weak suggestions of 
jet-cloud interactions in the optical emission line ratios. Typ- 
ically, we have found the emission line ratios in compact ra- 
dio sources to be remarkably similar to those in extended 
sources, with only a mild tendency for the optical line lu- 
minosities to be lower in compact radio sources for a given 
radio power (see Section 4.4). Possible explanations for these 
apparent contradictions include the following. 

• Difficulty in distinguishing between the models. 

As already noted in section 4.2 above, it can be difficult 
to distinguish in the diagnostic diagrams between the re- 
sults of simple slab AGN photoionization models and the 
shock+precursor models, and between the mixed medium 
photoionization models and the shock+precursor models, for 
some regions of parameter space. 

• AGN photionization masking shock signatures. 
It is entirely possible that the emission line clouds have been 
accelerated in jet-induced shocks, have cooled behind the 
shock, then been photoionized by the powerful AGN con- 
tinuum; the AGN may also photoionize the precursor gas. 
In this way the AGN may mask the ionizing effects of the 
jet-induced shocks. 

It will be necessary to obtain higher quality, spatially 
resolved emission line data (including temperature and den- 
sity diagnostics) for sources with well-characterised AGN, 
in order to truly distinguish the dominant ionization mech- 
anism for the warm gas in the compact radio sources. 

4.4 How do compact radio source compare to 
more extended radio sources? 

In addition to the data for this sample (see above), the di- 
agnostic diagrams in Figures [3] and [4] also show data for ex- 
tended radio radio sources, which includes nuclear regions 
and extended emission line regions, both with and without 
evidence for jet-cloud interactions (see Figure 3 for details) . 

As discussed in Section 4.2, the compact and extended 
radio sources occupy remarkably similar regions on the di- 
agnostic diagrams. This is a surprising result given the 
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Figure 5. continued. 

markedly different line profiles and more extreme kinemat- 
ics observed in the nuclear regions of compact radio sources 
compared to extended radio sources. There may be evidence 
in some of the diagrams (e.g. Figure 4b & 4c versus Fig- 
ures 3b & 3c) of a tendency for the broader, shifted compo- 
nents to be more consistent with the region occupied by the 
EELRs with evidence for jet-cloud interactions. However, 
our current sample of compact radio sources is too small to 
make any concrete statements. 

In order to compare the compact and extended radio 
sources further, we have also plotted radio -optical correla- 
tion p lots i n FigureO following, for example. iMorganti et al.l 
jl997h and lTadhunter et all (| 19981 ). Figure [5] includes all ob- 
jects from this sample, along with a redshi ft limited sub-set 
(0.05 < 2 < 0.7) of the 2Jy sample (e.g. iTadhunter et all 
Il99ct) . 

In Figure [S] the compact radio sources follow the same 
trends as the extended sources a nd occupy the high ra dio 
luminosity end of the correlations. iMorganti et alj l|l997h re- 
ported for their small sample (7 sources) that compact radio 
sources tended to have a lower [O IIIJA5007 luminosity at a 
given radio luminosity. This trend is also seen in our sample 
(Figure O, with the exception of two sources (3C 277.1 & 
3C 303.1) that lie toward the top boundary of the FR II 
trend. In addition, there is also evidence of this trend in the 
other lines plotted ([O IIJAA3727 and H/3). It is interesting 
that the compact radio sources tend to fall below the trend 
for extended sources, particularly as we would expect to ob- 
serve boosting in the optical flux if ionisation mechanisms 
other than AGN-photoionisation (e.g. jet-cloud interactions) 
were important. However, as discussed in Sections 3.3 & 4.1, 
the nuclear regions of compact radio sources are highly ex- 
tinguished and this may lead to the mild tendency of the 

3 Due to uncertainties in the reddening estimates (see Section 
3.3), the line fluxes used in this analysis are not corrected for 
reddening. 
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compact radio sources to lie below the trend observed for 
extended radio sources. 

Figure [5] also compares different emission line ratios 
with the radio luminosity. Although there appear to be 
no clear-cut differences between the compact and extended 
sources in terms of their line ratios, we note that the three 
GPS sources are amongst the group with the lowest [O II] /[O 
III] ratios. 



5 SUMMARY AND CONCLUSIONS 

It is clear that there are several trends regarding the physical 
conditions and ionisation mechanisms observed in compact 
radio sources. The main conclusions are listed below. 

• Dense and dusty circumnuclear cocoons. The ma- 
jority of the sources in this sample show evidence for high 
densities and large reddening in the nuclear regions, partic- 
ularly in the broader components. 7/14 sources show con- 
vincing evidence for increasing reddening with the line width 
of emission line component (for the remaining 7 sources, 4 
do not show the trend and we were unable to estimate the 
reddening for 3 sources), a st rong indication of a stratified 
ISM as in PKS 1345+12 fsee lHolt et ai1l2003h . This is also 
supported by an apparent trend in density with FWHM, al- 
though the uncertainties are larger. Hence, the results pre- 
sented here are entirely consistent with the idea of a dense 
and dusty cocoon of gas and dust enshrouding the nuclear 
regions, although evidence suggests that there is insufficient 
material to confine and frustrate the radio source. 

• Outflow driving mechanisms. The evidence for jet- 
driven outflows in the literature is strong. Extreme line 
broadening and outflow velocities are observed (up to ~ 2000 
km s _1 ) and there is strong alignment between the radio and 
optical line emission in all CSS sources with HST imaging. 
It is therefore surprising that our study has not provided the 
expected supporting evidence from the emission line ratios. 
Whilst our results suggest the gas densities and tempera- 
tures may be high, and hint at possible shock ionization in 
the broader, shifted components, the data do not clearly 
distinguish between shocks and AGN-photoionisation. 

• Quiescent, photoionised gas. As well as the kine- 
matically disturbed gas, we also observe quiescent gas on a 
variety of scales. In the nuclear regions, the narrow com- 
ponents of the emission line gas are predominantly pho- 
toionised (mainly mixed-medium models), consistent with 
observations of the nuclear regions of extended sources. In 
more extended regions we see evidence for both AGN pho- 
toionisation and also localised regions of star formation. 

• Similarities with extended sources. We present evi- 
dence for similarities between compact and extended sources 
suggesting that they are instrinsically similar objects. The 
high reddening in the nuclear regions is co nsistent with 
obser vations of, for example, Cygnus A (e.g. iTavlor et al.l 
2003). Also, the radio-optical correlations show that com- 
pact radio sources agree with the high radio power end of 
the correlations for extended sources, although there is a 
mild tendency for the compact radio sources to have lower 
optical emission line luminosities ([O II], [O III] & H/3) for a 
given radio luminosity. The scatter may be due to the high 
extinction in the nuclear regions of these sources. 



Hence, the evidence presented in this paper is fully con- 
sistent with the idea that compact radio sources form the 
early evolutionary stage in AGN with radio activity. 
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APPENDIX A: THE EMISSION LINE FLUXES 
(SAMPLE) 



Table Al: Emission line flux data for the sample. All line fluxes are quoted relative to H/3. The data are given for all apertures for each 
object, for all lines observed. Column 2 identifies the emission line component (n: narrow; i: intermediate; b: broad; vb: very broad) in 
accordance with the width definitions given in Section 2. For objects with a detected BLR, the data for this component is also presented 
and identified by BLR in the table. a Ha flux. b normalised to the Ha flux, f denotes lines obtained by 'free fitting' and are therefore 
the measured fluxes (1CF 16 erg s _1 cm~ 2 A -1 ) and not relative to H/3. 



3C 213.1 3C 268.3 


3C 268.3 


3C 268.3 


3C 277.1 


3C 277.1 


3C 277.1 


PKS 1345+12 


PKS 1345+12 


PKS 1345+12 


nuc nuc 


SE 


NW 


nuc 


SE 


SE HII 


nuc 


NW PA160 


SE PA160 



Mg II 2795.5 n 

i 

b 
vb 

[O III] 3133.70 n 80+10 
i 

b 
vb 

He I 3188.67 n 11+3 

i 

b 
vb 

[Nc V] 3345.8+ n 16+1 120+40 17+1 30+10 2+1 

20+20 

i 15+3 
b 
vb 

[Ne V] 3425.0+ n <20 <30 50+2 170+40 50+10 80+30 6+3 <286 

70+40 <82 

i 60 <20 50+20 

b 
vb 

[O II] 3727.64| n 850+30 190+10 110+4 650+40 160+10 350+50 200+10 205+5 1179+57 

220+10 170+90 102+16 

i 170 90+10 180+10 

b 
vb 

[Fe VII] 3759.99 n 13+2 
i 

b 
vb 



H/3 4860. 75+ 



0.9+0.1 



<0.8 



0.7+0.1 
0.9+0.1 
1.6+0.2 



0.97+0.03 0.05+0.02 



22+3 



8+3 



0.18+0.03 
0.13+0.02 



0.31+0.02 



2.58+0.14 



0.14+0.08 
0.49+0.07 



BLR: 100+10 
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APPENDIX B: APERTURE SELECTION AND 
EXTRACTED SPECTRA FOR ALL 
APERTURES (SAMPLE) 
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Figure Bl. Emission line spectra for 3C 213.1 along PA -61. 

Top: Spatial profile of the continuum and emission lines in 3C 213.1 along PA -80. The left panel shows the 2D spectrum of 
the most extended emission line along this PA, [O IIJAA3727 and the right panel shows spatial slices across the spectrum for both [O 
II]AA3727 (dotted line) and [O III] AA4959, 5007 (dashed line) emission lines. The fluxes in the different profiles have been scaled for 
comparison purposes. Overplotted are the positions and sizes of the apertures extracted. 

Bottom: Integrated rest-frame, continuum subtracted spectrum of the nuclear aperture. 



